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We discuss finite temperature phase diagrams of SU(N) gauge theory with massiess fermions as a 
function of the number of fermion flavors. Inside the conformal window we find a phase boundary 
separating two different conformal phases. Below the conformal window we find different phase 
structures depending on if the beta function of the theory has a first or higher order zero at the 
lower boundary of the conformal window. We also outline how the associated behaviors will help in 
distinguishing different types of theories using lattice simulations. 



There has recently been interest in studying the phase 
diagrams of massiess gauge theories, as a function of the 
number of colours, N, flavours Nf and fermion repre- 
sentations. This interest originates from applications of 
strong dynamics in elementary particle phenomenology, 
including unparticles [1] and (extended) technicolor sce- 
narios j2j-[4], but also from the purely theoretical aim of 
understanding the nonperturbative gauge theory dynam- 
ics from first principles. Recently several lattice collabo- 
rations have been pursuing this goal [5HT5] 



Several methods to estimate the vacuum phase dia- 
gram of a gauge theory exist. A traditional method is the 
ladder approximation to the Schwinger-Dyson equation 
for the fermion propagator, which yields an estimate for 
the critical coupling, a c — tt/(3C2(R), at the onset of chi- 
ral symmetry breaking. Comparing this with the value of 
the fixed point coupling a* obtained from solving for the 
zero of the two-loop beta-function /3(A) = — b X 2 — 61 A 3 , 
one can estimate the location of the boundary between 
conformal to confining phases as a function of N and Nf. 
[T71[TH]. To sketch a concrete example, let us assume an 
SU(iV) gauge theory with Nf fermions in the fundamen- 
tal representation, and define Xf = Nf/N c . In figure I] 
we show the vacuum phase diagram obtained under the 
ladder approximation as described above. 



as in QEEQ 

For applications in beyond Standard Model model 
building, one is mostly interested in the behaviors di- 
rectly below the conformal window; namely, how is con- 
formality lost in gauge theories. This has been discussed 
in [19] and more recently in [20]. Traditionally one as- 
sumes that there is a region where the theory exhibits 
quasi-conformality also known as walking. However, it 
may well be that such walking region does not exist. In 
essence, these behaviors are simple consequences of the 
properties of the zero of the beta function at the lower 
boundary of the conformal window. 
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The figure also includes schematic plots of the beta 
function of the theory in different portions of the phase 
diagram: for small values of Xf the gauge dynamics is 
QCD-like; the vacuum is confining and chiral symme- 
try breaking and the beta function is monotonic and 
negative. Within the conformal window, between the 
critical value x c ~ 4 and the value x as = 5.5, where 
asymptotic freedom is lost, the beta function has a non- 
trivial infrared fixed point (IRFP). Under the renor- 
malization group evolution, the coupling of these theo- 
ries shows asymptotic freedom at small distances, analo- 
gously to QCD, but flows to a fixed point in the infrared 
at large distances, where the theory hence looks confor- 
mal. Above the conformal window, Xf > x as the beta 
function is monotonic and positive and the coupling runs 



FIG. 1. An (N c ,Xf) phase diagram of SU(iV c ) gauge theory 
with Nf = N c Xf massiess fermions in the fundamental rep- 
resentation of the gauge group. The shaded region between 
the solid and dashed lines is the conformal window, while the 
lightly shaded region below the conformal window depicts the 
walking region which exists only if the zero of the beta func- 
tion at the lower boundary has a second (or higher) order 
zero. If the zero remains of first order all the way to the lower 
boundary there will not be walking (nor Miransky scaling). 

To concretize, concentrate first on a theory with walk- 
ing. To illustrate how the walking behavior arises, con- 



Of course, one cannot exclude the possibility of the existence of 
an UV stable fixed point here 
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sider the following form of the beta function, depending 
on an external parameter e controlling the approach to- 
wards the conformal limit: 

/3(A;e) = -c((A-A*) 2 -e), (1) 

where e > and the lower boundary of the conformal 
window is at e = 0. Then there are two scales, Air, at 
which the coupling diverges to infinity in the infrared, 
and Auv where the coupling runs towards zero. Between 
these two scales the coupling remains in the vicinity of the 
fixed point value, A ~ A*. Integrating the beta-function 
in ([lj, and applying these boundary conditions, one eas- 
ily derives^] 

^HX - cxp ( f XlR ~ (2) 

Am P lA- me))' 6 ' {2) 

so that the ultraviolet scale where the coupling runs to 
zero and the infrared scale where the coupling diverges to 
infinity become exponentially separated as e — > 0. Fur- 
thermore, all masses m oc Air, scale to zero exponentially 
when Xf approaches the conformal window from below. 

As the above calculation indicates, walking requires at 
least a second order zero in the beta function at Xf = x c . 
However, this is not the unique way how a generic gauge 
theory may exit or enter the conformal window. Consider 
a theory whose beta function is 

P(\) = -b (xf\ 2 -b 1 (xf\ 3 , (3) 

where &o ~ (a; as — xf remains positive until the value 
Xf — x^s i-e. until asymptotic freedom is lost at the upper 
boundary of the conformal window, while b± ~ (x c — xf 
changes sign at x c . This beta function is QCD-like all 
the way to x c when approaching the conformal window 
from below. Now, integrating this beta function we find 
that the infrared scale of the theory is 

A,»=<»(l + ^) , (4) 

where Ajj^ = no exp(— 1/(6qAo) is the limiting value as 
bi — > 0, i.e. at the lower boundary of the conformal win- 
dow. The central observation now is that the scale Air 
given above, decreases as one approaches the conformal 
window, but remains nonzero all the way to b = where 
it discontinuously drops to zero, i.e. its value inside the 
conformal window. 

A next logical direction is to consider these theories in 
finite temperature. Like vacuum phase diagrams, under- 
standing the finite temperature phase diagrams of generic 



2 One can also consider more precise beta function including also 
QCD-like behavior at A -> 0, like /3(A) = -cA 2 ((A- A*) 2 -e); the 
main result, i.e. the Miransky scaling, remains as in the simple 
example we consider here. 



gauge theories will help us to understand the dynamics of 
strong interactions. The phenomenological applications 
on the other hand include phase transitions in early uni- 
verse (13 [52]. So far less analytic efforts, and hardly 
any lattice efforts, have been devoted to expose the finite 
temperature phase diagrams. In this paper we point out 
that there are several interesting features to be uncov- 
ered; our central result is the identification of qualitative 
structures shown in Figs. [2]and[3j The figures show the 
possible phases of SU(N) gauge theory as a function of 
temperature and Xf = Nf/N c . Fig. [5] corresponds to 
theories which exhibit walking behavior while Fig. [3] cor- 
responds to theories where there is no walking. In other 
words, these figure correspond, respectively, to theories 
where the zero of the beta function at the lower bound- 
ary of the conformal window is of second or first order. 
In both cases, for small Xf, the finite temperature phases 
are analogous to a QCD-like theory, while inside the con- 
formal window (i.e. Xf > x c ) there is a phase transition 
between two different conformal phases. The order of the 
finite temperature phase transitions in each of these fig- 
ures depend on the underlying theory. Across all phase 
boundaries, the order parameter is provided by the free 
energy itself; i.e. there will be a (substantial) change in 
the number of effective degrees of freedom. We will now 
discuss different regions of the phase diagram in detail. 

The difference between the two cases shown in Figs. 
[2] and [3] arises when approaching the conformal window 
from below. First, if the theory exhibits walking for val- 
ues x* < Xf < x c , this is reflected at finite temperature 
as an existence of a quasi conformal intermediate phase 
where p/T A ~ const, over a temperature range whose 
width increases rapidly as the conformal window is ap- 
proached. At Xf = x c and T — there is a second order 
quantum phase transition. Directly below the conformal 
window at some fixed Xf such, that x* < Xf < x c the 
system is in the confined hadronic phase at low temper- 
atures. Then, as the system is heated, there is a phase 
transition (or a smooth crossover) to the quasi conformal 
phase at T^jr ~ Air, which extends to T" c .uv ~ Auv- 
The scales Air uv correspond with the ones defined in 
([2]), and quasi conformality arises due to the fact that 
/3(A) ~ 0, i.e A ~ const, between A IR and Auv- The 
masses of the lightest hadronic states are nih ~ 27tT C) ir, 
and as x — > x c these vanish due to Miransky scaling. 
Above T c .uv the system is in the asymptotic partonic 
plasma phase. We note that these behaviors have been 
observed also in an holographic model for walking gauge 
theory thermodynamics [231 121] • However, if the theory 
does not have walking behavior, the phase diagram cor- 
responds with the one shown in Fig. [3j In this case the 
zero temperature transition at Xf — x c is of first order, 
and the finite temperature behavior below the conformal 
window is QCD-like all the way to x — ¥ x c . There is a 
single deconfining phase transition at T c ~ Air, where 
the IR scale is as defined in Q. 



3 



The robust features, present in both cases, are: First, 
at small values of Xf the theory is expected to display 
behavior akin to QCD; confinement at low temperatures 
and a high temperature plasma of collectively free par- 
tons. These behaviors are well understood. For funda- 
mental representation matter there is a single deconfin- 
ing phase transition likely driven by restoration of chiral 
symmetry [26] . With higher representations the possibil- 
ities are richer since deconfinement and chiral restoration 
are not necessarily strongly intertwined 



SU(N) gauge theory, massless fermions 
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FIG. 2. An (x, T) phase diagram of gauge theory with mass- 
less fermions where the dynamics of the theory provides walk- 
ing behavior. 



SU(N) gauge theory, massless fermions 
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FIG. 3. An (x, T) phase diagram of gauge theory with mass- 
less fermions but without walking behavior. 

Second, inside the conformal window x > x c there is 
a transition between two conformal phases: from a low 
temperature IR phase of "unparticles" to deconfined par- 
tonic plasma at high temperatures. To show how this 
behavior arises, consider first the Banks-Zaks region, i.e. 
values of Xf close to upper boundary of the conformal 



window, say Xf — x as - e, e < 1. Then the fixed point is 
perturbative [31], and we expect perturbation theory to 
be applicable for the whole evolution between the ultra- 
violet and infrared scales. To be explicit, we consider 
SU(A r ) gauge theory with fundamental fermions, and 
identify A = Ng 2 /(8ir 2 ) where g 2 is the gauge theory 
coupling constant. Pcrturbatively, the pressure is given 
by [321133] 



p(T) = (ai + a 2 87r 2 A(T)/A0T 4 



(5) 
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We apply these formulas with the two-loop coupling A(/i) 
determined by numerically solving for its evolution from 



/3(A) 



dX 
din fi 2 



= -b a X 2 - 6xA 3 



(6) 



Here the first two coefficients of the beta function are 
b = 11/3 - 4/37/a;/ and h = 34/6 - 2C f /NT f x f - 
20 /6TfXf. For fundamental representation the required 
group theory factors are T/ = 1/2 and C/ = (./V 2 — 
l)/(2iV). In computing the pressure we choose the renor- 
malization scale /i = 2nT. The resulting phase diagram 
is shown in Fig. [4] for SU(3) gauge theory with 16 funda- 
mental representation fermions. In this case the transi- 
tion is not a sharp phase transition, but a continuous 
cross over which takes place over a wide temperature 
range. This is due to the fact that the magnitude of 
the coupling and the beta-function is very small over the 
entire range of couplings, < A < A* < 1. One can still 
assign a " critical temperature" , defined as a maximum 
of the interaction measure (c — 3p)/T 4 ; i.e. the trace of 
the energy momentum tensor. The results corresponding 
the the pressures in Fig. [4] are shown in Fig [5] 

These results imply that, as Xf is lowered from the up- 
per boundary of the conformal window, T c (xf) decreases, 
and the difference between the effective degrees of free- 
dom in the low temperature and high temperature phases 
increases. 

What can we then say about the finite temperature be- 
haviors towards the lower boundary of the conformal win- 
dow? Here the properties of the low temperature phase 
are dictated by strong coupling physics, but some robust 
results can still be derived due to conformal invariance. 
Consider the trace anomaly of the underlying gauge the- 
ory. In the massless limit it is 



(7) 



in terms of the coupling and the beta-function of the 
theory. On general grounds we expect that for conformal 
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where A = A(/io) Taking the thermal expectation value, 
i@u) = 6 ~~ 3p> choosing the scale as fx = T and applying 
standard thermodynamic relations, we obtain 
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FIG. 4. Pressure in perturbation theory for SU(3) gauge 
theory. Different curves correspond, from top down, to 
Nf = 16, 15, 14, 13 flavors in the fundamental representation. 
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FIG. 5. The interaction measure, (e — 3p)/T , in SU(3) gauge 
theory with Nf = 16, 15, 14, 13 fermions in the fundamen- 
tal representation. The critical temperature, identified with 
the maximum of the interaction measure decreases as Nf de- 
creases. 
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where 7 is the anomalous dimension. The nontrivial tem- 
perature dependence then arises as the anomalous dimen- 
sion changes as the theory is heated and evolves under 
the renormalization group between the IR and UV fixed 
points. To determine the leading temperature behavior 
of (8^) in the deep infrared, we linearize the beta-function 
near the IRFP at A = A*, 



(9(A) ~ c(A - A*), c>0. 
This can be easily integrated to yield 



A( M ) = A* + (A - A*) 



Mo 



Here C is a constant of integration and reflects the IR 
T/GeV degrees of freedom, C ~ jir. We therefore predict a sim- 
ilar phase structure as in the perturbative Banks-Zaks 
region: At very low temperatures p/T A is a constant and 
increases as T is increased. In high temperatures we ex- 
pect p/T 4 to saturate to the Stefan-Boltzmann limit cor- 
responding to the hot gas of free quarks and gluons. In 
high temperatures the running of the coupling is gov- 
erned by the UVFP at A = 0, and the temperature de- 
pendence can be described by perturbative results. One 
can hence derive limiting behaviors in the deep infrared 
and ultraviolet domains, but first principle methods are 
needed to bridge the two. 

Having established the qualitative features of finite 
temperature phase diagrams of gauge theories far, near or 
inside the conformal window, we now discuss the conse- 
quences in light of future lattice simulations. For theories 
outside the conformal window, i.e. confining gauge theo- 
ries, an obvious goal would be to determine the value of 
x* where the quasi conformal phase should appear and 
phase structure changes from QCD-like theory with a 
single finite temperature phase transition. In particular, 
there would be Miransky scaling exhibited by the critical 
temperature associated with the deconfining transition 
from the hadronic phase. However, mapping out the ex- 
tent of the quasi-conformal intermediate phase at larger- 
values of Xf might be numerically more challenging prob- 
lem due to large scale separation which imposes the need 
for large lattices. Furthermore, it is important to realize 
that it may well be that the theory does not have walk- 
ing behavior at all, and in such a case there would not 
appear the quasi conformal phase either. 

Therefore, it might be interesting to start inside the 
conformal window, where it should be possible to deter- 
mine the phase boundary between the low temperature 
unparticle phase and high temperature partonic plasma 
phase. The lattice results, in combination with analytic 
insights obtained so far should help in drawing a coherent 
picture. For example there is currently a consensus that 
the simple SU(2) gauge theory with two Dirac fermions 
in the adjoint representation of SU(2) has an IRFP con- 
trolling the vacuum properties [101 E] . Heating up 
this theory would provide an excellent testbed for estab- 
lishing the existence of the two distinct conformal phases 
at finite temperature, as we have outlined in this work. 
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